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Abstract 



Vortex-line equilibrium configurations and critical currents for type-II su- 
perconducting films at zero-temperature are studied theoretically. The films 
are assumed to be of thickness less than or equal to the penetration depth, free 
of pinning by imperfections, and to be subjected to a magnetic field parallel 
to the film surfaces and to a transport current perpendicular to the field. By 
numerical minimization of the exact London theory energy expression, using 
simulated annealing techniques, the equilibrium configurations are determined 
with great accuracy over a wide range of fields and currents. These consist of 
chains of straight vortex-lines whose number depends both on the field and 
on the current. Transitions involving a change by one in the chain number 
are found to take place either if the field is increased at zero transport current 
or if the transport current is increased at constant field. At the transitions, 
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there is a considerable rearrangement in the vortex-line positions and a small, 

discontinuous, change in their number. The equilibrium configurations with 

chain numbers that are not too small form a nearly perfect triangular lattice, 

centered with respect to the film surfaces. However, small deviations from 

this arrangement are found to be important in determining the behavior near 

the transition and when the transport current approaches the critical current. 

It is found that the critical current has a a non-monotonic dependence on the 

field. The zero-temperature equilibrium phase diagram in the field-transport 

current plane is reported. 
74.76-w; 74.60-w 
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I. INTRODUCTION 



There is a great deal of activity nowadays in the investigation of the properties of the 
mixed-state in type-II superconductors under various physical situations. One problem of 
interest is the study of vortices in films of thickness comparable with the penetration depth 
placed on an external magnetic field parallel to the film surfaces. In this case vortices enter 
the film as straight lines parallel to the field direction. The shielding current penetrate 
the full thickness of the film and pins the vortex-lines to the film's interior. A transport 
current applied parallel to the film surfaces and perpendicular to the vortex-lines will not 
dissipate energy unless its magnitude is large enough for the Lorentz force exerted by it on the 
vortex-lines to overcome the force generated by the shielding currents and by the interactions 
between the vortex-lines. Therefore, even in the absence of pinning by imperfections, the 
mixed state in the film can support a finite transport current without dissipating energy, 
that is in thermodynamic equilibrium. 

Previous theoretical work on this problem has calculated, in the London limit, the equi- 
librium vortex-line configurations and critical currents using both analytical |]-[5[ and nu- 
merical |j) methods . These find that the equilibrium configurations consist of vortex-line 
chains and that there are transitions, induced by changing the applied magnetic field, in 
which the number of chains changes. 

Measurements on films in this geometry of the magnetization and of the critical current 
have been reported by several authors ||6Hl2|. It is found that the magnetization has peaks 
as a function of the field, not present in the bulk material, which are interpreted as resulting 
from the transitions involving the change in the number of chains . Some of the critical 
current measurements find non-monotonic dependence on the field that is also attributed to 
these transitions P-pH • Other measurements find no such dependence, but the large critical 
currents observed are interpreted as due to pinning of the vortex-lines by the shielding 
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The aim of this paper is to calculate in detail the vortex-lines equilibrium configurations 
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and critical currents at zero temperature in films of thickness less than or equal to the 
penetration depth, without pinning by imperfections, in the presence of both a parallel field 
and a transport current. The starting point of this calculation is the known expression for 
the energy of a system of straight vortex-lines parallel to the film surfaces, exact in the 
London limit for sufficiently thin films Minimization of this energy is carried out 

numerically, using simulated annealing techniques, assuming only that the distributions of 
vortex lines in the film is periodic. By this method the equilibrium vortex-line configurations 
are determined very accurately over a wide range of external fields and transport currents. 
Several new results are reported. 

The equilibrium configurations obtained by this method are found to consist of vortex- 
lines chains, all with the same intra-chain spacing between vortex-lines. In the absence of 
a transport current, it is found that, as the external field grows above the lower critical 
field, there are regions where the number of chains remains constant in which the number of 
vortex-lines in the film grows smoothly with the field. Transitions consisting in a change by 
one of the number of chains are found to take place at critical field values. These results are 
consistent with previous ones, obtained using a similar numerical method |]]. The present 
calculation reveals that at the transitions, besides the considerable rearrangement in the 
vortex-lines positions found previously ||, their number has a small discontinuous increase. 
Even for a modest number of chains, the equilibrium vortex-line configurations are found to 
be close to a triangular lattice, as predicted in Ref. . However, it is found that there 
are significant differences between the triangular lattice and the equilibrium configurations, 
not reported previously, that play an important role in the vicinity of the above described 
transitions. The changes in the equilibrium vortex-line configurations caused by a transport 
current are studied in detail and the critical current is calculated. These results are also new. 
It is found that, by changing the transport current at constant field up to the critical current, 
two distinct behaviors take place, depending on the field value. In the first the number of 
chains remains the same as that for no current. The effects of the current are to shift the 
chains towards the film surface and to increase a little the number of vortex-lines. Near 
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the critical current this shift is non-uniform, contrary to previous calculations that find an 
uniform shift [Q]. It is found that, if the vortex chains were forced to shift uniformly, a critical 
current considerable larger than that obtained from the exact calculation would result. The 
second type of behavior is one in which, as the current grows past a certain value, smaller 
than the critical current, a transition to a configuration with one additional chain takes 
place. This is similar to the transition that occurs by increasing the field without a transport 
current and has not been predicted in previous calculations. The critical current versus field 
curve is found to be non-monotonic with a structure that is a direct consequence of the 
vortex-line chains configurations that exist in the film. The zero-temperature equilibrium 
phase diagram in the external field-transport current plane is obtained. 

This paper is organized as follows. In Sec. ||. the London theory results for the energy 
are briefly reviewed. In Sec. |T| the numerical method for energy minimization is presented. 
In Sec. [TV] the results obtained by this method are reported and their physical significance 
is discussed. Finally, in Sec. [V] the predictions of this paper are compared to experiment 
and its conclusions stated. 

II. LONDON THEORY 

Consider a film made of a type-II superconducting material of thickness D, length L 
and width L. The material is assumed to be uniaxially anisotropic, with anisotropy axis 
perpendicular to the film surfaces (the c-axis) and characterized by anisotropy parameter 
7 (= m a b I m c ) . The other material parameters are the penetration depth A, for currents 
flowing parallel to the a-b plane, and the coherence lengths in the a-b plane £. 

The film is assumed to be subjected to an external field H parallel to the film surfaces 
(along the b-direction) and to a transport current flowing perpendicular to H and also 
parallel to the film surfaces (a-direction), whose average current density is Jt ( Fig. [I]). 

The interest here is to obtain the equilibrium vortex-line configurations and the critical 
currents at zero temperature for D < A, in the London limit, A > (. It is assumed that: 
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i) the vortex-lines are straight and parallel to the field direction (b-direction). ii) Their 
distribution in the ac-plane is periodic in the a— direction, with period a v . This periodic 
structure consists of L/a v identical primitive unit cells with dimensions Dxa v placed side by 
side. Each cell contains n v vortex-lines located at (x a j, x C j), j = 1, 2, n v , with < x a j < a v 
and < x C j < D ( Fig. [I]). Even though these assumptions allow for very general vortex-line 
configurations, the calculations reported here find that the equilibrium ones consist of n v 
chains of vortex-lines located at {x C j}, with vortex lines uniformly spaced from one another 
within each chain by a v ( Fig. |l|). 

The solutions of London equations for the film with an arbitrary distribution of straight 
vortex-lines parallel to the film surface were obtained by several authors using the method 
of images One such a vortex-line, with vorticity q = 1 and located at (x a ,x c ) (D > 

x c > 0), generates two infinite sets of mirror images. One set has images with vorticity —q 
located at — x c + 2nD, n = 0, ±1, ±2, .... The other has images with vorticity q located at 
x c + 2mD, m = ±1,±2, .... 

The total energy of an arbitrary distribution of such vortex-lines can be written as 

E = E int + E sel f + E H + Ej . (1) 

In Eq. ([I]) E int is the energy of interaction of vortex-lines with one another and with the 
images, excluding the self-images; E se if is the vortex-lines self energy plus the energy of 
interaction with their own images; Eh is the energy of interaction of the vortex-lines with 
the screening current generated by the external field and Ej is the energy of interaction of 
the vortex-lines with the transport current. 

For an isotropic film (7 = 1), and in the limit where (nX/D) 2 ^> 1, E int and E se if are 
given by 

2 1 ^ cosh 7r (na v + x ai - x aj )/D - cos 7r(x ci + x cj )/D 

E mt /L 2 = e— ^ In — ■ J -f— J -f—] , (2) 

2a " nXj cosn K(na v + x ai - x aj )/D - cos n(x ci - x cj )/D 

where e = ((/>o/47rA) 2 and J2 n ,i,j runs over n = 0, ±1, ±L/a„, and i,j = l,2,...,n v , 
excluding for n = the i = j term ; 
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Esds/L = h ln [ ww ] ' {) 

The energies of interaction with the currents are given by 

Air\ 2 H ^ cosh(x CJ -D/2)/A 

EH/L=e ^7^} coshD/2A 1] ' (4) 

and 

2 8n 2 \ 2 DJ t sinh(x c ,-D/2)/A 

Ej/L " 6 o a,c £ [ sinhD/2A + 1] ' (5) 

The total energy for an anisotropic film (7 / 1) is related to that for an isotropic film 
(7 = 1), with the same A and D, by the following scaling relation |L3| 



E[£, a v , {(x aj , x cj )}, H, J t ; 7] = jE^, ja v , {(>yx ajl x cj )}, 7 if, 7 J t ; 7 = 1] . (6) 

Thus the equilibrium configurations for an anisotropic superconductor follow from the cor- 
responding ones for the isotropic material by scaling the parameters according to Eq. (^j). 
In what follows the discussions will be limited to isotropic materials. 

In order to obtain the equilibrium vortex-line configuration it is necessary to find the 
minimum of E, Eq. ([I]), for given H and Jt- This is discused next. 

III. ENERGY MINIMIZATION 

The energy minimization is carried out numerically using simulated annealing. For a 
given H and Jt = 0, the following procedure is used. Assuming that there are n v vortex- 
lines within the unit cell, a standard Monte Carlo (MC) procedure that moves vortex-lines 
within a single cell is used to find the positions (x a j,x C j) and cell length a v at a given 



"temperature" [|14]]. This fictitious temperature is used to overcome energy barriers. By 
lowering it the minimum of E/L 2 , Eq. (|IJ), is found. The equilibrium configurations are 
obtained from this data by comparing the minimum energies for different n v and determining 
the lowest one. 
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For J t 7^ it is necessary to search for a local minimum of E/L 2 , Eq. (|l|), at fixed n v . 
The reason is that, for fixed n v and for J t > 0, E is unbound because vortex-lines placed at 
x c = D have zero interaction energy and self-energy, but negative Ej. Thus the energy can 
always be lowered by putting more and more vortex-lines at x c = D close together. The 
local minimum is, clearly, the physical relevant one. 

For given H and n v the local minimum is determined as follows. Starting from the 
vortex-line configuration found for these values of H, n v and for J t = 0, a new configu- 
ration that minimize E is determined for small J t by running the MC procedure at very 
low "temperature". This avoids "thermal" fluctuations that might overcome the energy 
barrier separating the local minimum from the unphysical configurations described above. 
By further increasing J t by small amounts new local minima configurations are found. The 
equilibrium configurations for J t > are determined as in the case where J t = 0, described 
above, by comparing the energies of the local minima for different n v . It is found that above 
a certain value of J t the vortex-lines accumulate at x c = D, their number growing without 
limit. This mimics what happens above the critical current, where vortex-lines enter the 
film at x c = 0, are dragged by the transport current to the film surface at x c = D, where 
they are annihilated. The smallest value of of J t where this happens is interpreted as the 
critical current J c . 

This method allows the equilibrium configurations of interest to be determined accurately 
using small values of n v and with modest personal computers. 

In the results reported in Sec. |V| A and £ are considered as fixed lengths. Films with 
£ = 10~ 2 A and D = A, A/2, A/4 are investigated over a wide range of values of H and J t . 
According to Eq. (|]), for these values of D, E is weakly dependent on A. The mathematical 
derivation of Eq. (|2]) requires only that (ttX/D) 2 ^> 1, which is satisfied. 

The natural units for the external field and for the transport current are the lower critical 
field, H c i, and of the depairing current, Jd, respectively. For the films with D < A these 
quantities are given by 
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0o ln(2£)/7rfl 
cl A 2 4tt(1 - l/cosh£>/2A) ' 1 J 

and 

T - 1 m 
d ~ 12^ £A 2 ' K ) 

For the films studied here the values of H& in units of (po/X 2 are: 
H cl (D = A) = 2.92, H c i(D = A/2) = 9.04, and F c i(.D = A/4) = 28.37. 



IV. RESULTS AND DISCUSSION 

For J t = 0, and for H just above H&, the equilibrium configuration is a single vortex- 
chain (n v = 1) located at the film center ( Figs. § and |3|). As H increases the vortex-lines 
in the chain come close together, that is a v decreases, as shown in Fig. |j. Eventually a value 
of H is reached where the single chain splits in two chains (n v = 2), located symmetrically 
about the film center ( Figs. ^| and |3|). This two-chain configuration exists over a range of 
H values, above which the energy minimum corresponds to three chains (n v = 3) ( Figs. 
and [3].). As H is further increased the pattern is as follows. A stable n^-chain configuration 
occurs over a range of H, followed by a n v + 1-chain configuration as H grows past a critical 
value. At the n v —>■ n v + 1 -transition a v increases discontinuously, with the amplitude of 
the discontinuity decreasing with increasing n v ( or H), as shown in Fig. |]. Similar results 
are obtained for other D values. 

In Fig. H the total number of vortex-lines per unit length of the film n v /a v is shown as 
a function of H . It is seen that the number of vortex-lines increases roughly linearly with 
H, with constant slope in the regions where n v is constant. At the n v — > n v + 1 chains 
transition there is a discontinuous jump in n v /a v , by about 10% in the region covered by 
Fig. [5|. The relative size of the discontinuities in both a v and n v are considerably larger ( 
Fig. |J). However, these quantities increase by about the same factor so that the relative 
increase in their ratio is smaller. At higher values of H/H c \ than those shown in Fig. |5| the 
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number of vortex-lines continues to increase in the same fashion, but the discontinuity at 
the transitions become smaller. 

Even for modest values of n v , the equilibrium configurations are found to be very close to 



of Eq. (HD show important differences with this triangular lattice, as discussed next. 

In the centered triangular lattice the distance between adjacent chains is d c = D/(n v +l), 
which is also the distance from the first and last chains to the film surfaces at x c = and 
x c = D, respectively, and a v = d c j sin 60 ( Fig. ||). The results of the present calculation show 
that the chain spacing increases slightly towards the film center. Numerical minimization 
of E, Eq. (H), assuming a centered triangular lattice structure, gives the a v x H curve with 
steps shown in Fig. |j. In this approximation a v ~ D/(n v sin 60), and n v depends on H 
through the minimization of E/L 2 . Consequently, a v is independent of H in the range 
where n v is constant, in sharp disagreement with the results described above. However, it is 
apparent from Fig. [| that the triangular lattice prediction agrees well with the mean value 
of a v in the range where n v is constant. At the n v — > n v + 1 chains transition d c and a v 
decrease discontinuously to fit a centered triangular lattice with n v + 1 chains. This leads 
to a (n v /a v ) x H curve with steps: n v /a v is constant in the range where both n v and a v are 
constant and increases discontinuously by A(n v /a v ) ~ 2n 1 ,sin60/D at the transition. This 
linear dependence on n v of A(n v /a v ) is also in strong disagreement with the results of the 
full minimization of Eq. ([!]), discussed above. 

A simpler way to estimate n v and a v for large follows. First replace the dis- 

crete vortex distribution by a continuous one z/,;(x c )(=number of vortex-lines per unit area). 
Minimization of Eq. ([I]) under this assumption gives 



This result shows that the vortex-line density is slightly larger near the film's surfaces, in 
accordance with the results of the full minimization of Eq. ([!]), as discussed above. Next 



a perfect triangular lattice centered with respect to the film's surfaces ( Fig. |]), as predicted 
in Refs. . However, the vortex line configurations obtained by numerical minimization 




(9) 



10 



approximate the average vortex-line density, n v /a v D, by the average over the film thickness 
of u c (x s ) which, according to Eq. ([)]) is z/ c = (H/<f) )(2X/D) tanh (D/2X), and use for a v 
the centered triangular lattice approximation discussed above. For n v 3> 1 this gives n v = 
D/(a v sin 60) and 

r 2XH sin 60 , D. 1/2 
a ' = t^0^ tanh 2A> ' ■ (10) 

This prediction agrees well with the average of a v in the interval where n v is constant, even 
for modest values of n v , as show in Fig. It shows that this average value decreases with 
H^ 1 / 2 . However, this approximation overestimates n v in the range shown in Fig. |4]. 
For D < X the magnetic induction 

0o ^ h cosh (x cj -D/2)/X 1 
B = ^D^ coshD/2A ] ' (U) 



resulting from the vortex-line configurations in the field is small. It follows from Eq. ( |TT|) 
that B ~ <p (n v /a v D)(D/X) 2 / 12. Consequently, M = (B — H)/4ir deviate very little from 
the Meissner state value M = —H/kn. 

For J t > 0, and if H is not too close to the field where the n v — > n v + 1 chains transition 
takes place at Jt = 0, the equilibrium vortex configuration has the same n v as that for 
J t = and the same H. The chains positions are shifted towards the film surface located 
at x c = D, a v decreases slightly and, consequently, the number of vortex-lines in the film 
increases a little. The shift grows with J t and is non-uniform, with the chain closest to film 
surface at x c = D shifting the most, as shown in Fig. [7] . In this figure it is seen that as 
Jt Jc the chain closest to the x c = D film surface is about halfway between its position 
at J t = and this surface. 

A simple model to calculate J c , proposed in Ref. , is to assume that the vortex-lines 
form a rigid triangular lattice. In this case the effect of the force exerted by the transport 
current on the vortex-lines is to displace the lattice uniformly towards the film surface. 
Equilibrium is reached when this force is balanced by that exerted by the other vortices 
and by the shielding current. Numerical calculations based on this model and on Eq. ([[]) 
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show that above a critical value of Jt the rigid lattice has no equilibrium position within the 
film. However, this critical value is found to be considerably larger than that calculated as 
previously described. This shows that, at least for the range of H values studied here, the 
subtle changes in the vortex-lines configuration away from the triangular lattice caused by 
J t play an important role in determining the critical current. 

If H is close to the field for the n v — > n v + 1 chains transition at Jt = 0, it is found that 
a similar transition takes place by increasing J t at constant H. 

The vortex-line structure in the film can be summ arized in a zero-temperature equilib- 
rium phase diagram, as shown in Fig. ||. For D = A/2 and D = A/4 the phase diagram is 
similar. The critical currents for these D values are show as a function of H in Fig. |9]. 

V. CONCLUSIONS 

The calculations reported in Sec. [IV] describe the behavior of an ideal model. That is, 
vortex-lines in equilibrium at zero temperature in films without imperfections under a field 
that is precisely parallel to the film surfaces. Even a very small angle between the field 
and the film surfaces, which is always present in experiments, invalidates these predictions, 
since dissipation occurs when J t ^ 0. In such fields the vortex-lines are tilted. The Lorentz 
force exerted by J t on these lines has a component along the b-direction which will move 
the lines and dissipate energy. However, since the tilt is small, it takes only a small amount 
of pinning to prevent this motion. Morever, it is reasonable to assume that the tilted 
vortex-line configurations are very close to those for untilted vortex-lines at the same field. 
That is, the vortex-lines are nearly parallel to the b-direction and their average positions in 
the ac-plane coincides with those obtained for lines parallel to the b-direction. Thus, the 
ideal model predictions are applicable to equilibrium properties of films, provided only that 
pinning by imperfections is sufficiently weak and that the field makes a small angle with the 
film surfaces. 

It has been usual to relate peaks in the magnetization data on films with the n v — > n v + 1 
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chains transitions that occur in the ideal model |6|-|8[. It is therefore interesting to discuss 
what are the predictions of the ideal model to the behavior of the magnetization. 

As discussed in Sec. [TV], the magnetization due to vortex-lines that are exactly parallel 
to the film surfaces deviates very little from the Meissner value. 

If the vortex-lines are tilted with respect to the b-direction the equilibrium magnetic 
induction has a c-component B c . Assuming that the tilt is uniform, B c is proportional 
to the number of vortex lines per unit area, that is B c oc <f)o(n v /a v D). The equilibrium 
magnetization has also a c-component M c oc <fio(n v /a v D — H/<fio). Thus, M c is a direct 
measure of number of vortex-lines in the film. In Fig. |5| (inset) m c = —(n v /a v D — H/<fro)\ 2 
is plotted versus H. The peaks seen in this figure result from the increase in the number of 
vortex-lines in the film that take place at the 1 — > 2 and 2 — > 3 chains transitions. Similar 
results are obtained for larger values of D. Experiments on Nb/Cu multilayers reported 
in Refs. find that the M c x H curve has peaks at field values approximately equal 

to those where the n v —>■ n v + 1 chains transition takes place in the ideal model. If these 
experiments measure the equilibrium M c , there is the following disagreement with the ideal 
model predictions. The experimental | M c | decreases in the range of H where n v is constant, 
whereas the ideal model predicts that it increases. If in the experimental setup the vortex- 
lines in the film are not in equilibrium, then it is unclear how the M c peaks are related to a 
transition that only takes place in equilibrium. 

In the experiments reported in Refs. M c is measured as H is changed at J t = 0. 

According to the results reported above the n v — > n v + 1 chains transition can also take place 
if Jt is increased at constant H . This suggests that, if peaks in M c are indeed related to the 
n v — > n v + 1 chains transitions, it may be possible to observer peaks in M c by increasing J t 
at constant H. 

The measurements of the critical current in films under a parallel field reported in Refs. 
P^T2| can also be compared to the ideal model predictions. The main characteristics of these 
predictions are a non-monotonic J c x H dependence and values of J c ~ 1 — 9 x Jd- The 
results reported in Refs. |p|-[TT|] show non-monotonic dependence on H, with a characteristic 
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features at H values comparable to those where n v — > n v + 1-chains transition take place. 
However the critical currents measured in Refs. P,p"Tf are one order of magnitude smaller 



than those predicted by the present calculation. The results reported in Ref. show only 
a monotonic dependence on H, and find J c values that are about one order of magnitude 
larger than those shown in Figs. [5] and 

In conclusion then, the ideal model predictions are not fully supported by the exper- 
imental data. However, this model is useful as a starting point to study the interesting 
behavior of films with weak pinning under parallel fields. Further work is necessary in order 
to incorporate into the model the effects that are neglected here. 
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FIGURES 

FIG. 1. A portion of the film cross-section showing a typical vortex- line configuration. The 
primitive unit cell is indicated by the double arrow. The field, H, is entering the paper ( 
b-direction). The transport current direction is indicated by the single arrow. 

FIG. 2. Equilibrium configurations for J t = and D = A. The field values are: H = 2H c i for 
n v = 1; H = 3H cl for n v = 2; H = 5H cl for n v = 3; H = 7.25H cl for n v = 4; H = llH cl for 
n v = 5 and H = l5H c i for n v = 6. 

FIG. 3. Vortex-chain positions for J t = and D = A. The continuous lines are just guides to 
the eye. 

FIG. 4. Intra-chain vortex-line spacing for Jt = and D = A. Thick line: full minimization of 
E/L 2 . Dashed curve: continuous approximation. Inset: comparison with the centered triangular 
lattice approximation (thin line). Both approximations are described in the text. 

FIG. 5. Number of vortex-lines per unit length in the film for Jt = 0. Inset: magnetization 
along the c-direction for D = A/4 and Jt = 0. 

FIG. 6. Circles: vortex-line configuration for H = 30H c i, n v = 9, Jt = and D = A, ob- 
tained from full minimization of E/L 2 . Triangles: centered triangular lattice for n v = 9. The 
configurations are symmetric with respect to x c /D = 0.5 

FIG. 7. Vortex-chain positions for D = A and H = 28H C \ ( n v = 9) . 

FIG. 8. Zero-temperature equilibrium phase diagram for D = A. 

FIG. 9. Critical currents. The arrows on the upper (lower) axis mark the range over which n v 
has the value indicated for D = A/4 (D = A/2). 
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